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ABSTRACT 

The scale length over which convection mixes mass in a star can be calculated as the inverse 
of the vertical derivative of the unidirectional (up or down) mass flux. This is related to the 
mixing length in the mixing length theory of stellar convection. We give the ratio of mass 
mixing length to pressure scale height for a grid of 3D surface convection simulations, covering 
from 4300 K to 6900 K on the main-sequence, and up to giants at logg = 2.2, all for solar 
composition. These simulations also confirm what is already known from solar simulations, 
that convection doesn't proceed by discrete convective elements, but rather as a continuous, 
slow, smooth, warm upflow and turbulent, entropy deficient, fast down drafts. This convective 
topology also results in mixing on a scale as that of the classic mixing length formulation, and 
is simply a consequence of mass conservation on flows in a stratified atmosphere. 

Subject headings: Stars: convection — mixing length — methods: numerical 

1. Introduction 



For stellar structure and evolution calculations (e.g., tbenl 119671) . models for the relevant energy- 
transport mechanisms are needed. In the interior, radiative or conductive transfer of energy is trivial, al- 
though the atomic physics going into the opacities and conductivities are not. A satisfactory description 
of convection in a manner amenable to stellar structu re calculations, however, has prov en elusive. Stellar 



modelers usually use the mixing length theory (MLT) ( B o hm- Vitensell 1 95 a : iGoughll 1 9771) . which despite its 



simplicity and many shortcomings, has been highly successful in matching a variety of stellar observations. 
In this model, convection is thought of as a collection of fluid parcels that travel a distance equal to the mix- 
ing length and then thoroughly mix with the surrounding ambient medium, sharing their energy. The MLT 
picture is completely symmetric, so the upward flow of buoyant warm elements is accompanied by an equal 
number of cooler and over-dense elements moving inward. In these calculations the mixing length is a free 
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parameter and is specified as a multiple of the pressure scale height Hp = Pj pg, where P is the pressure, 
p is the density and g is the gravitational acceleration. In general, there are also free parameters specifying 
the geometry of the fluid parcel - the ratio of its horizontal to vertical size and another specifying the lateral 
radiative transfer between convective elements and the ambient medium. The simple mixing length model 
is used due to a lack of better alternatives and because it is quick to compute. Its primary function is to 
determine the entropy of the deep, efficient and nearly adiabatic part of the convection zone in the star. If the 
magnitude of the mixing length and the other parameters could be specified a priori, simple mixing length 
theory would be capable of fulfilling this function. However, there is little guidance from theory as to their 
values. 

In addition, the mixing length model suffers from several other problems: it does not allow for over- 
shooting motions that mix material into the surrounding stable layers, it does not allow for turbulence, it 
does not allow for asymmetry between upflows and downflows, especially the large temperature difference 
near the surface and its effect on opacities and radiative transfer. These are all phenomena that, according to 
realistic 3D simulations, contribute to the atmospheric entropy jump and hence, the adiabat of the deep, ef- 
ficient convection. Several attempts have been made to impro ve the mixing length m odel. Some have taken 



the mixing length to be the distance to the nearest boundary (IStothers & Chin 



1997), some have introduced 



a two-stream model with separate properties of the upflows and downflows dNordlundl 119761) . some have 
introduced the idea of a spectrum of sizes of the convective cells (ICanuto & Mazziteililll99ll) . 



Several attempts have been made to determine the appropriate value of the ratio of mixing length to 
pressure scale height by matching theoretical and observed evolutionary tracks in the Hertzsprung-Russell 



diagram (Chieffi et al 



1995 



Stassun et al 



2004h. bv matchin g observed and model propert ies of binary 



system s ( Fernandes et al. 19981 20021 : IMontalban et al. 2004 ). and red giant temperatures dFerraro et al. 



2006). There have also been a few attempts to determine the mixing length parameter by comparing 
one-dimensional mixing length stellar model calculations with two- and three-dimensional numerical sim- 
ulations of stellar convection by either matc hing entropy in the adiabatic portion of the convection zone 
dLudwigetal.lll999l: ^reytaa & Salaris 199 9|) or by matching t emperature and density at a given pressure in 



the iriterior ^raffl oedach et al. 



1999 



Trampedach et al.lbOllbh . 



Mixing length theory works reasonably well for determining the thermal structure of convection zones 
because real mixing of mass occurs in stratified convection. Stars are stratified. Their density decreases from 
their center to their surface in order to produce a pressure gradient to balance the inward pull of gravity. In 
order to conserve mass in such a stratified med ium, most of the ascendi ng material must turnover and head 
back down within about a density scale height (IStein & Nordlundlll989l) . Thus, mixing really does occur on 
the order of a density scale height. 

It is possible to easily and unambiguously calculate the mass mixing length from computational fluid 
dynamic simulations of stellar convection. Such simulations, when they employ realistic equations of state 
to relate the pressure to the density and internal energy or temperature and employ sufficiently accurate treat- 
ment of radiative transfer to determine the heating and cooling, have as free parameters only the resolution 
of the computational grid used and the magnitude of the diffusion used to stabilize the calculations. These 
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parameters, however, are not free to be adjusted to match observations. The diffusion is minimized against 
the requirement of a numerically stable simulation, and the resolution is chosen as a trade-off between re- 
solving features in the photosphere and computational cost. We would not lower the resolution or increase 
the diffusion to give a better match to some observations. 

In this paper, we present results of such mass mixing calculation for a homogeneous sample of convec- 
tion simulations, covering a large part of the Hertzsprung-Russell diagram for which stars have convective 
envelopes. The simulations form an irregular grid, covering 4 500-6 900 K on the main sequence and going 
up to giants with log (7 = 2.2 and T e g =4 300-5 00 K, all with solar compos ition. The simulati ons em- 



but with many improvements. The simulations are described in more detail in lTrampedach et all (|201 laUb|) . 



ploy the Mihalas-Hummer-Dappen equation of state ( Hummer & Mihalasl ll988). lin e-opacity from Kurucz 
(1992) and continuum opacities largely from the Marcs stellar atmosphere pac kage (Gustafsson et al.lll975h . 



2. Mass Mixing Length 




Pressure (dyn /cm 2 ) 

Fig. 1. — Uni-directional downward mass flux as a function of the mean pressure for the stars in our sample. 
The break in the mass flux near pressures of 10 5 dyn/cm 2 occurs in the photosphere. The effect of the bottom 
boundary can be seen as a small downturn near the high pressure end of the traces. 

The mass mixing length, the distance over which, because of mass conservation, most of the ascending 
fluid turns over and starts to descend, is the inverse of the logarithmic derivative of the unidirectional mass 
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flux dStein et alJl2009h . 



IdlnF^^/drr 1 
Idlni^CrJ/drr 1 . (1) 



This quantity is readily determined from three-dimensional, stellar convection simulations. The equality of 
the mass flux in the up- and downflows, is just another way of stating mass conservation at each depth. 

In the rest of this paper we will look at the upflow, which mass flux is 

F tm (r) = (p(r)V(r)) t A t (r) « (p(r)h(V(r)) t A t (r) , (2) 

The brackets, (f(r))f, denote horizontal averages over the upflows of a quantity /: (/(r))f = f(s, r)d 2 s/A^(r), 

with A^(r) being the area with upflows and s is the coordinates in the horizontal directions. The approxi- 
mation of the right-hand-side of Eq. © ignores correlations between density and velocity, but we only use 
this approximation to gain some insights into the behavior of £ m (r), by re- writing it as 



din (/>)-)- + dln(F) t + dlnvl t 



dr dr dr 



(3) 



Below the surface region, the area of the upflows is nearly constant at about 2/3 of the total area, so its 
gradient has little effect on £ rn . From Eq. (0 we see that, ignoring gradients in velocity and filling factor, 
the mass mixing length is simply the density scale-height, H p = dr/d In p. Another way to look at it, is 
to consider an upflow along the decreasing density of an average stratification with gradient If the 

upflow had retained its density, it would be overdense by Ap/p = — ^^Ar, compared to its surroundings, 
after having traveled the distance Ar (remember that the gradient is negative). This is obviously unstable, 
and instead the fraction Ap/p overturns into the downflows, in order to keep the hydrostatic balance. That 
fraction is 1 for Ar = 1/^ = H p in a linear stratification and 1/e in an exponential stratification. Locally 
the upflow will therefore be "eroded" with an e-folding scale of H p . This "erosion" of the upflows does not 
mean that they loose coherency and disappear over that length-scale — on the contrary. We observe coherent 
upflows over the entire 5.3-7.3 pressure scale heights covered by the convection zones of our simulations (a 
similar number of pressure scale heights are included above the convection zones of the simulations). We 
likewise see downdrafts that reach all the way to the bottom of our simulations (in all of our simulations), 
with weaker ones either stopping or merging into fewer strong downdrafts, before they reach the bottom. 
The mass mixing length simply describes how only a very small fraction of the mass in the upflows make 
it from the interior and all the way up to the photosphere, simply due to mass conservation in a stratified 
atmosphere. 

Realistically the velocity decreases inward from the top of the convection zone, so its derivative is of 
opposite sign to that of the density. This decreases the variation of the unidirectional mass flux and increases 
its scale length compared to the density scale height. In a deep solar simulation, we find a logarithmic 
derivative of the upward mass flux of ~ 4/(5H p ) (Figure [3]), and the vertical velocity must therefore vary 
approximately as 

Kcrt OC p" 1 / 5 . (4) 
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The various terms of Eq. © are shown for the deep solar simulation in Figure [2] We clearly see that 
the density gradient is the main term and that the velocity gradient produces the deviation from l m = 
H p . The gradient of the filling factor is only significant in the surface layers. The difference between the 
approximation of the right-hand-side of Eq. d3) and the full expression, can be seen by comparison with 
Figure |3j which is based on the latter. 
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Fig. 2. — The three terms of Eq. © constituting the mass mixing lengths, in terms of the density scale 
height, Hp (left), and the pressure scale height, Hp (right). The full expression is shown with solid line. 
The effect of ignoring the gradient in filling factor is shown with dot-dashed, and just including the density 
gradient results in the dashed line. The example is shown for the deep solar simulation, also presented in 

Rg. m 



Near the surface, the mass mixing length is not a constant multiple of either the pressure or density scale 
heights (Figs. [2] and [3]). This is due to the abrupt changes in all three mass-flux ingredients in the surface 
layers. The density gradient becomes very small (for all but the coolest dwarfs the density in the upflows has 
a sub-photospheric inversion) resulting in a spike in the density's contribution to the mixing length (See right 
hand panel of Fig. 0). The small density gradient is however accompanied by a large velocity gradient due 
to buoyancy braking of the upflows where the density gradient becomes small or inverted, which reduces the 
density induced peak in the mixing length. In the interior the density and velocity gradients have opposite 
signs and a nearly constant ratio. However, as the surface is approached the buoyancy driving increases, 
which steepens the velocity gradient and gives rise to the wide sup-photospheric bump in the mixing length. 
The photospheric side of this bump is considerably smoothed by the contribution from the filling factor. 
Division by the pressure scale height to get the l/Hp presented in Figs.|4]|5]decreases the difference between 
the atmospheric and the interior values. 

To determine the mass mixing length, i m , for a simulation, we first calculate the horizontal and time 
averages of the vertical upward (and downward) momentum (which is the unidirectional mass flux). The 
unidirectional mass flux increases with increasing depth as the density increases and the velocity decreases, 
but more slowly than p~ 1 / 2 . We then calculate the inverse of the logarithmic derivative of this total up- 



-6- 



or downward unidirectional mass flux. We have computed th e mass mixing length for the set of realis- 
tic, compressible, convection simulations by iTrampedach et all (1201 lbl) . Figure Q] shows the unidirectional 
(downward) mass flux for all the stars in our sample. Next we determine the multiple of the pressure and 
density scale heights that best match the inverse logarithmic derivative of the unidirectional mass flux be- 
tween its local maximum at the steep temperature gradient near the surface and where the bottom boundary 
conditions begin to affect it, due to boundary effects on the velocity, seen as the upturn at the bottom in 
Figs. [3]-[5] This increase in the mass mixing length near the bottom boundary is an artifact of the bottom 
boundary condition that the upflows are forced towards vertical at the boundary. 



3. Results 



l/H ? 



.70±0.09 
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log 10 (P/P e //) 
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log, (^e»f) 
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log 10 (P/P„ f ) 



Fig. 3. — Sun's mass mixing length (solid) in units of the pressure scale height (left), density scale height 
(middle) and distance to the top of the convection zone (right) as a function of the mean pressure in units of 
the surface pressure. Dashed line is the average and dotted lines are the extremes away from the boundaries. 



For the Sun we have deep (20 Mm) convection calculations (IStein et al. 2009) and the determination of 
the mass mixing length is straightforward, as shown in Figure [3] From this we see that the ratio of the mass 
mixing length to the pressure scale height, ap, is nearly constant with depth in the interior of the convection 
zone until the bottom boundary is approached. The ratio to the density scale height, a p , is closer to one, but 
varies slightly and linearly with depth. The mixing length is greater than the density scale height because 
buoyancy is continually accelerating the flow resulting in a significant contribution from the second term in 
Eq. ([3]). The mass mixing length in units of the distance to the top of the convection zone, a r - rcz , is also 
nearly constant with depth in the interior, but has a value of only 0.45 rather than one. In all three cases the 
constant proportionality is destroyed at the same point near the bottom boundary due to boundary effects 
(See Fig. [T]). At the top of the convection zone, in the highly superadiabatic layer, the mass mixing length 
becomes large compared to all three measures. Hence, larger values of the mixing length parameters, ax, 
are needed near the surface than in the interior, wi th a r _ rc7 becoming largest of all, a factor of 5, compared 
to only a factor of 1.5 for the two other methods. ICanuto & Mazzitellil (119911) suggest that a mixing length 
equal to the distance to the boundary, av-r cz = 1, should work for the entire Sun. This does not match the 
actual mass mixing in the simulations which is half of the expected value in the interior, adiabatic part of the 
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convection zone, and a factor of five larger near the surface. From a relative logarithmic pressure of 0.5 and 
down to 4.7 where the bottom boundary effects are felt, this normalization of the mass mixing length results 
in a concave curve, and it is not obvious which range should be considered close to constant. If we allow 
larger fluctuations around the mean, that mean value will increase slightly as will the range over which the 
ratio is within these fluctuations. To obtain a range similar to that of ap would require accepting fluctuations 
as large as ±0.27, apart from them not being random fluctuations, but rather a systematic concaveness. 

The reasoning behind choosing £ = r — rcz is based on the MLT concept of convective elements having 
size £, and that they cannot be larger than their distance from the top of the convection zone. Convection 
simulations show, however, that the convective topology is one of upward and downward moving streams 
with continual entrainment of upward moving fluid into the downflows. Convection is a continuous not a 
discrete process. There is therefore no conceptual problem with an t > r — rcz- 

As ap is the most constant measure of the mass mixing length in the interior of the deep solar simula- 
tion (Figure [3]>, and we have similar experiences with other simulations, we advocate a mixing length for- 
mulation using ap. We note that with an adiabatic stratification in hydrostatic equilibrium, Hp = H p /ji. 
So Hp < Hp due to the simultaneous variation of p and T, and consequently ap ~ 71 a p . In the deep 
solar simulation 71 reaches the completely ionized value of 5/3, but for the other 26 simulations, listed in 
Table [Q only the six simulations along the high-T c fj/low-g edge get close to fully ionized, and the rest do 
not get closer than 71 ~1.22-1.32 (their minima are in the 1.16-1.22 range). The radial variation of 71 
therefore also affects the differences between mixing length normalizations with H p or Hp, although the 
small gradient in the interior makes this a minor effect. 

For other stars the available simulations are shallower (about 3 dex less in pressure), so some of the 
simulations have a smaller range between the surface effects and the bottom boundary effects, and the 
mixing length for those simulations are, for the moment, less precisely determined. 

Table Q] presents the mass mixing length in units of Hp for the 26 stars for which we have performed 
convection simulations and for which the simulations have a deep enough domain to determine the mass 
mixing length. The table is ordered in order of increasing log g, and for similar gravity, in order of increasing 
T c g. Figure @]-[5] show the individual ap determinations and Figure [6] shows how ap varies as a function of 
effective temperature and gravity. 

4. Discussion and Summary 

The basic structure of stellar convection zones is determined by the need to conserve mass. Stars are 
highly stratified so most of the mass ascending from depth must turn over and go back down within the 
order of a scale height. The high density mass from the interior can not be supported in the low density 
outer layers. This produces a true mixing by entrainment of the almost laminar upflows into the turbulent 
downflows. There is therefore a real physical basis for the mixing length theory of convection. From 
computational fluid dynamic simulations of stellar convection zones one can calculate this real mixing of 
mass to determine the length scale on which it occurs. Naively, one might expect it to be some multiple of 



-8- 



Table 1 : Mass mixing length and fundamental parameters for the 26 simulations 



Simulation 


JV1K class 


T cff /[K] 


log g 


71 A~ 1 71 jT 

M/Mq 


l/Hp 


1 


K3 


A f f A 

4669 


O O AA 

2.200 


1 A CA 

3.450 


"I f 1 A 1 A A r" 

1.69 ± 0.05 


2 


K2 


4962 


1 O A A 

2.200 


A 1 AA 

4.300 


i tt i a r\r* 

1.77 ± 0.06 


3 


1/ C 

K5 


4294 


ni a o a 

2.420 


A O A O 

0.842 


1 1 A A T 

1.75 ± 0.07 


4 


Kl 


4994 


2.930 


2.350 


1.71 ± 0.07 


5 


G8 


C c r/ 

5556 


3.000 


r» CAA 

2.500 


"1 f T 1 A A O 

1.67 ± 0.08 


6 


i r / \ 

K0 


CO o o 

5288 


1 /I O 1 

3.421 


1.833 


1 TO 1 A AO 

1.73 ± 0.03 


-7 

7 


K3 


4718 


1 CAA 

3.500 


A A A A 

0.940 


-i ry f | A A f 

1.75 ± 0.05 


o 

8 


K0 


pi OA 

5189 


1 CAA 

3.500 


1.660 


"1 TO 1 A AO 

1.72 ± 0.03 


9 


F9 


6111 


3.500 


1.748 


1 TA 1 A AT 

1.70 ± 0.07 


1 A 

10 


G6 


5671 


1 A/11 

3.943 


1.110 


"1 ^7f7 1 A a r - 

1.77 ± 0.05 


1 1 


Procyon F4 


/^c o a 

6529 


3.966 


1 CAA 

1.500 


"1 Z" 1 T 1 A A f 

1.67 ± 0.06 


1 ni 

12 


K3 


4673 


/I AAA 

4.000 


0.727 


"1 AA 1 A AO 

1.90 ± 0.08 


13 


K2 


4980 


4.000 


0.825 


1.85 ± 0.15 


14 


F4 


6611 


4.000 


1.490 


1.68 ± 0.06 


15 


F9 


6147 


/I A /I A 

4.040 


1 1 OA 

1.180 


-1 TO 1 A AA^ 1 

1.72 ± 0.006 


16 


Gl 


cao a 

5929 


4.295 


1 AO O 

1.028 


"1 TA 1 A AO 

1.70 ± 0.02 


1 "7 

17 


K4 


/I Al 

4603 


/I O AA 

4.300 


A C AO 

0.502 


"1 (\f 1 A A/ -1 

1.96 ± 0.06 


18 


K0 


5327 


4.300 


0.788 


1.83 ±0.05 


19 


F2 


6908 


4.300 


1.396 


1.72 ±0.05 


20 


Sun G5 


5778 


4.438 


1.000 


1.76 ±0.08 


21 


K4 


4500 


4.500 


0.553 


2.05 ± 0.05 


22 


F7 


6286 


4.500 


1.201 


1.69 ±0.03 


23 


Kl 


5021 


4.550 


0.767 


1.98 ±0.26 


24 


G2 


5894 


4.550 


1.085 


1.76 ±0.04 


25 


G9 


5480 


4.557 


0.933 


1.85 ±0.04 


26 


K4 


4531 


4.740 


0.800 


2.20 ± 0.04 



the density scale height. However, as we observed for the Sun in Figure [3] it is actually the ratio of the mass 
mixing length to the pressure scale height which is nearly constant over many orders of magnitude change 
in the pressure. All of our values of mass mixing lengths are therefore presented in terms of the pressure 
scale height, Hp. 

In Figure [6] we present the variation of a = IjHp with the atmospheric parameters T e g and log g. 
First off, comparing the range of a with the uncertainties in Table [Q we see that the variation over the HR 
diagram is small, but signific ant. The range of values found here (1.6-2.2) is consistent with those used in 



stellar evolution calcu l ations (Salaris & Cassisi 2008 ; Montalba n~& D'Antonall2006 : Montalban et al. hood : 



Fernandes et alJI 19981 : IChieffi et al. 



1995 



Pedersen et al.lll990h . although such calculations assume a con- 



stant value for a as the star evolves. We have also compared the mass mixing length values found here with 
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0.0 0.5 1.0 1.5 2.0 2.5 0.0 0.5 1.0 1.5 2.0 2.5 0.0 0.5 1.0 1.5 2.0 2.5 

log 10 (P/Peff) lO gl0 (P/P e „) log 10 (P/P ef f) 

Fig. 4. — Ratio of mass mixing length to pressure scale height for stars labeled by their simulation number 
in the table. Solid line is the calculated mass mixing length divided by the pressure scale height. Dashed 
lines are the minimum and maximum range chosen for the ratio a = I /Hp. 



those determined by matching the adiabat i n the convection zone fr om the simulations with that obtained 
from mixing length theory stellar models (ITrampedach et al.ll2011bh . We find excellent agreement except 
for the coolest, low mass, dwarf stars where the mass mixing length values are larger. It should also be re- 
membered that in addition to the other mixing length parameters specifying geometry and radiative energy 
exchange, other input physics, atmospheric opacities and the T-r relation, also influence the appropriate 
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4 



13: a P =1.85±0.15 



14: a P =1.68±0.06 




15: a p =1.72±0.006 



5 
4 



0.0 0.5 1.0 1.5 2.0 2.5 0.0 0.5 1.0 1.5 2.0 2.5 

lOg 10 (P/P e ff) lOg 10 (P/P e ff) 



Fig. 5. — Next set of stars (See Fig.@]for explanation). 



6000 5000 
WW 



Fig. 6. — The mass mixing-length in units of Hp f or the 26 simulations , as listed in Table [Jin the T e g - log 



g plane. We have also plotted evolutionary tracks (ISchaller et al.lll992uCharbonnel et al. 
as indicated at the bottom of each track. 



1999) for masses 



value of the mixing length to use in stellar structure and evolution calculations. 

We find a is smallest towards the high-T c ff and low-g side of the diagram. These are also the sim- 
ulations that have the least efficient convection, where high speeds and large superadiabatic gradients are 
necessary for carrying the flux to the photosphere. The least efficient convection has the smallest mass 
mixing lengths, as expected. The most efficient mixing occurs in the coolest dwarf with the most sedate 
convection. The solar value from our grid of simulations of a = 1.76 ± 0.08 agrees well with several 
determinations from the literature, lending credence to the solar models and to the common practice of 
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calibrating «mlt against such models. We also see from Figure [6] however, that the solar value is less ap- 
plicable for lower mass stars, which also have the largest change of a during their evolution. The evolution 
of higher-mass stars and the ascent up the Hayashi-track, seems well described by a single value of a. 

The surface layers cannot be properly accounted for with a MLT-like formulation. Crucial compo- 
nents like over-shooting, turbulent pressure, radiative transfer in the presence of large temperature inho- 
mogeneities, etc., are not included in this simple model, but greatly affect the surface layers. The mass 
mixing length found here, is the appropriate mixing length to use in a formulation that can realistically in- 
clude all these phenomena — an ideal MLT formulation. Still lacking such a formulation, the solar mass 
mixing length found in the present work, does not necessarily agree with that found from matching solar 
evolution models to the Sun's present luminosity and radius. The latter furthermore depends on all the other 
ingredients going into the atmosphere of the evolution model, such as atmospheric opacities and how the 
photospheric transition from optically thick to optically thin is treated, e.g., through some match to an at- 
mosphere model or through T-r relations. The mixing length of MLT is strongly connected to the actual 
length-scale of mixing, as found in this paper, but this connection is weakened by all these confounding 
factors entering into the MLT version. 

The increase in the ratio of mass mixing length to pressure scale height in the superadiabatic layers 
near the surface raises the question of why stellar evolution calculations with the interior values of a work 
well in matching observations of global stellar properties. In this regard, note that to match line profiles (e.g. 
Balmer lines ) a small value of a that produces a steep te mperature gradient is ne eded in the atmosphere. In 



this case, a large value of a must be used in the interior (|Montalban et al.ll2004|) . Our result suggest such a 



two value set of mixing length parameters may indeed be appropriate. However, the main conclusion should 
be that mixing length theory is not well defined and different variations are possible that will yield similar 
results. 

Calculations, such as presented here, will reduce the uncertainty in the appropriate choice of mixing 
length in stellar structure and evolution calculations. In particular, we have used the behavior with depth, 
to discriminate between various mixing length formulatio ns, and find that the mass mixing length, I, is best 



described with a constant ap = £ /Hp. In a future paper ( Trampedach et al.ll201 lbh we present calibrations 



of ap(MLT) against the simulations, based on ID MLT envelope models that employ the exact same EOS 
and opacities as the simulations, and u se a new self-consistent treatment of T-t relations, computed from the 
simulations ( Trampedach et aliboila ). That calibration will be aimed at reproducing the asymptotic adiabat 



with a standard MLT formulation, and will contain our recommended values for ap (MLT). Our present cal- 
culation of a is model independent in that we compute the actual physical mixing length of the simulations, 
without connecting it to a particular MLT formulation. Our present result, that a p seems convergent with 
depth, is a confirmation that MLT is a physically reasonable model for stellar convection, despite its many 
shortcomings, par ticularly near the surface . The present work also confirms the soundness of the scheme 



employed in the iTrampedach et all (|201 lbh ap (MLT) calibration, since, for most of the simulations, it re- 
sults in a-values very similar to what we find here. This agreement between the two independent methods, 
suggests that ap(MLT) has been approximately disentangled from the various uncertainties in atmospheric 
physics that have historically been absorbed into «mlt- 
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With the results from iTrampedach et al.l (1201 laUbh . supported by our present results, ID MLT models 
can be made more self-consistent and with a reduced set of free parameters. From an MLT point of view, 
there are of course still the undetermined geometric form factors that could conceivably be adjusted to match 
the adiabat in conjunction with the mass mixing length found here. 
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NASA grant NNX08AI57G. RFS was supported by NASA grants NNX07AO71G, NNX07AH79G and 
NNX08AH44G, and NSF grant AST0605738. This support is greatly appreciated. 
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